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FIGURE CAPTIONS 


Figure 1. Model of a Strike-Slip Fault in a Layer Lying Over a Half-Space. 

Figure 2a, Rheological Element for Standard Viscoelastic Solid. 

2b. Rlieological Element for Maxwell Substance, 

Figure 3. Displacements Parallel to the Fault, u(t), for Model DILE; = I m; Tick Mark Spacing = 2L 
Figure 4. Displacements Perpendicular to the Fault, v(t), for Model DILE; u^ = 1 m. 

Figure 5. Shear Strain, Mfodel DILE; ii^ = 1 m. 

Figure 6. Normal Strain, eji (t), for Model DILE; u^ = 1 m. 

Figure 7. Normal Strain, for Model DILE;Uq - ! m. 

Figure 8. Rotation co(t) for Models DILE and DILV; u^ « I m. 

Figure 9. Displacements Parallel to the Fault, u(t), for Model DILV; Uq = 1 in. 

Figure 10. Displacements Perpendicular to the Fault, v(t), for Model DILV; Ujj = I m. 

Figure 1 L Shear Strain, Gjjit), for Model DILV; Uy = 1 m. 

Figure 12. Normal Strain, for Model DILV; Uy = 1 m. 

Figure 13. Normal Strain, £32 Model DILV; Uy = 1 m. 

Figure !4a. Parallel Displacement, u(t), Versus Distance from the Faiilt for Model DILE; Uy = 1 m, 

L = 50 kin for tlie Scales on Top and Right Side of Figure, 

1 4b. Parallel Displacement u(t), Versus Distance from the Fault for Infinite Fault; Uy = 1 m, 

H = SO km for tlie Scales on Top and Right Side of Figure. 

Figure 15a. Shear Strain, ej 2 (t), Versus Distance from the Fault for Model DILE; Uy = 1 m, L = 

50 km for the Scales on Top and Right Side of Figure. 

!5b. Shear Strain, e, 2 (t), Versus Distance from the Fault for Infinite Fault; Uy = 1 m, L = 

50 km for the Scales on Top and Riglit Side of Figure. 
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ABSTRACT 

This paper, the first of two related articles, presents a model of 
viscoelastic deformations associated with earthquakes. A strike-slip fault 
is represented by a rectangular dislocation in a viscoelastic layer (litho- 
sphere) lying over a viscoelastic half-space (asthenospherc). Deformations 
occur on three time scales. The initial response is governed by the instan- 
taneous elastic properties of the earth. A slower response is associated 
with viscoelastic relaxation of the lithosphere and a yet slower response is 
due to viscoelastic relaxation of the asthenospliere. The major conceptual 
contribution of this paper is the inclusion of lithospheric viscoelastic pro- 
perties into a dislocation model of earthquake related deformations and 
stresses. Numerical calculations using typical fault parameters reveal that 
the postseismic displacements and strains are small compared to the co- 
seismic ones near the fault, but become significant further away. More- 
over, the directional sense of the deformations attributable to the elastic 
response, the lithospheric viscoelastic softening, and the asthenosphcric 
viscoelastic flow may differ and depend on location and model details. 

The results and theoretical arguiricnts suggest that the stress changes 
accompanying lithospheric relaxation may also be in a different sense 
than and be larger than tlie strain changes. 


POSTSEISMIC VISCOELASTIC SURFACE DEFORMATION AND STRESS 
Part 1 : Theoretical Considerations, Displacement and Strain Calculations 


INTRODUCTION 

The subject of tliis study is the spatial and temporal dependence of postseismic displacements, 
strains, and stresses. In this paper we develop a theoretical model of the postseismic deformations 
and report on calculations of surface displacements and strains. In a companion paper (Cohen and 
Cook, to be published), herein referred to as Part 2, we report the corresponding stresses and exa- 
mine the effects of variation in fault parameters on the computed deformations. 

Our analysis is based on a viscoelastic model of the earth’s rheology and employs features from 
dislocation theory (Steketee, 1958a, 1958b; Chinnery, 1961 , 1963;Press, 1965; Rybicki, 1971) 
and the correspondence principle (Flugge, 1967). The assumed model for a strike-slip fault is shown 
in Figure 1 . The lithosphere is modeled as a layer of thickness H lying over an asthenospheric half- 
space. The strike-slip faulting occurs uniformly over a rectangular surface from -L to L and from 
some upper level which is not necessarily the surface to a depth D witliin the lithosphere. As will be 
seen below, the deformations depend on both the deviatoric and dilatational rheology of the litho- 
sphere and the deviatoric rheology of the asthenosphere. We assume that both the lithosphere and 
asthenosphere exhibit viscoelastic deviatoric responses although the details of the response and the 
associated time scales are different. Tlie dilatational response of the lithosphere is taken to be elastic 
in one version of our model and viscoelastic in another. Specifically we model the lithospheric devi- 
atoric behavior with the standard linear solid shown in Figure 2a. The dilatational response is repre- 
sented by an elastic element in model DILE and a standard linear solid in model DILV. In the latter 
case the deviatoric and dilatational time constants are assumed equal. The deviatoric response of the 
asthenosphere is modeled by a Maxwell clement shown in Figure 2b. 


T]\e time constants associated with the iitliospheric and astlienospheric viscoelasticity are ex- 
pected to be markedly different. Although there is some debate about the physical mechanism for 
the lithospheric viscoelasticity we take an empirical approach based on laboratory and field investi- 
gations which yield time constants ranging from several minutes to a few months (e.g. Robertson, 
1964). On the other hand astlienospheric viscoelasticity has a longer time constant expected to 
range from years to decades (e.g. Nur and Mavko, 1974), or longer. 

Two more preliminary comments need to be made about this model. First the consideration 
of the lithosphere’s dilatational response is important only for finite length faults. For infinitely 
long faults only the deviatoric behavior has to be considered and the model is two-dimensional. We 
have already reported some results from this model (Cohen, 1979) and will discuss others below. 

The present results comprise the more detailed study to which we alluded in the earlier paper. Sec- 
ond, as the thickness of the lithosphere increases (with fault depth held constant), the influence of 
the asthenosphere decreases. In the limit that the lithosphere becomes infinitely thick, model DILV 
gives no post-seismic defonnations but model DILE does. This general result, based on the corres- 
pondence principle, is discussed in many references. Viscoelastic calculations based on an infinite 
half-space model of the earth have been reported by Singh and Rosenman (1974), and Rosenman 
and Singh (1973a, 1973b). Models involving dislocation theory and a layer over a half space have 
been studied by Rundle and Jackson (1977), Savage and Prescott (1978), and Rundle (1978). A 
viscoelastic model of the Palmdale bulge has recently been reported by Rundle and Thatcher (1979). 

THEORETICAL CONSIDERATIONS 

The starting point for our analysis is the set of equations describing the static displacements 
associated with a dislocation in an elastic layer lying over an elastic half-space. A set of approxi- 
mate equations have been derived using an image technique by Rundle and Jackson (1977) who find 
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where the symbols used in these equations are defined on Table I and where (Chinnery, 1961, 
1963) 
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In order to find tlie corresponding time dependent viscoelastic displacements we can follow 
the procedure dictated by the correspondence principle. We write the stress-strain equations for 
the viscoelastic material as 


Pa = Qe 



where P and Q are linear operators. For the' case of pure elasticity, P = 1 and Q = 3k or Q = 2/i 
depending on whether dilatational or dcviatoric behavior is under study. In the more general case 
wlien P and Q are linear operators we take the Laplace transform of equation 4 to find 

P(s) o = QCs) T (5) 

where P and Q are polynomials in the Laplace transform variabJe s. Equatioji 5 shows that the 
Laplace transformed viscoelastic solution can be obtained from the elastic solution by making 
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It is also necessary to include the transformation 

— U- 

Up = -Si since Uq = a constant for t > 0 (7) 

The resulting equations for u and v can then be subject to an inverse transformation to deter- 
mine the viscoelastic displacements. This procedure is very cumbersome. A more expeditious 
procedure results from two simplifications. First we impose the physically reasonable condition 
that the lithosphere relaxes rapidly compared to the asthenosphere. Next we restrict the discus- 
sion to the three times of greatest interest, namely, t = tg = 0 — the time immediately after the 
earthquake, t = tjg — a time long compared to that required for lithospheric relaxation, and t = 
ta — a time long compared to that required for asthenospheric relaxation. In these three time 
limits the displacement equations become elastic with effective reduced moduli , P 2 > and k as 
shown on Table II. The parameter substitutions shown on Table II are obtained directly from con- 
sideration of the viscoelastic stress-strain equations. Consider, for example, the deviatoric equation 
for a three element standard linear solid 
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For a rapid change in stress and strain (e.g. during an earthquake) the terms in o and i dominate 
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Thus for rapid changes 


The two viscoelastic models that we have considered have in common the standard linear solid re- 
presentation of the dcviatoric behavior of the lithosphere and the Maxwell deviatoric behavior 
of the asthenosphere. They differ in that model DILE takes the lithosphere to be elastic in 
dilatation and model DILV takes to be a standard solid with Pa/Pb = and rilii^ = r?7kb* 
The displacements at the earth’s surface (z = 0) are in these models; 
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where the coefficients aj and bj depend on time and the model; they are shown on Table III. The 
time dependent quantity V has been evaluated in an earlier paper (Cohen, 1979) and is given by 
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Tlie strains are obtained from appropriate differentiations on the displacements; 
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From these we obtain tlie shear strain 
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The elastic displacements and strains depend on the location (x, y), the source parameters 
(AU,D, L, H) and the ratio of elastic constants k (or k^) to and to p. It is convenient to 
choose as a unit of spatial measurement the half-fault length, L, and to normalize the results to a 


unit dislocation, AU. Furthermore as is common in theoretical geophysics we take the instanta- 


neous) Poisson’s ratio to be corresponding to = — and set Pi(to) = P’(to). This 

^ii(to) 3 


leaves as parameters the quantities D/L and H/L and the coordinates x/L and y/L, The choice 
of equal instantaneous shear moduli for the lithosphere and asthenosphere sets F(tQ) = 0 and all 
terms under the summation signs in equations 14, 15, and 17-20 vanish in the calculation of the 
initial elastic deformations. The deformations associated with the lithospheric relaxations depend 
on the parameter /3 of the relaxed standard solid and the parameter P (tg) which compares the 
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effective moduli of tlie lithosphere and asthenospliere. Tlie deformations occuring as the astheno- 
sphere relaxes are independent of the asthenospheric modulus as this half-space has no long-term 
rigidity. 

We will defer until part 2 a detailed discussion of the stresses that accompany these defor- 
mations. However it is worth noting that the surface stresses may have significantly different 
time liistories from the strains. For example, the shear stress, a^ 2 , is given by 
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As time progresses from t^ to tg not only does the strain change, but so does the effective shear 
modulus coupling the stress to the strain. Since the fractional change in ei 2 is small in many 
cases and since 0 may have a value of order unity, the modulus change can dominate the early 
postseisniic stress history. In fact the reduction In the effective shear modulus causes a decrease 
in the initial shear stress drop associated with the eartiiquake. As a consequence there is a 
partial stress recovery at some points even thougli the strain undergoes a postseismic decrease. 

A similar effect has been seen sliding block models of earthquakes that use standard linear solid 
viscoelastic elements (Dieterich, 1972; Cohen, 1978). 

RESULTS 

In this section we will discuss some figures wliich illustrate the deformations predicted by 
the equations presented above. For this discussion we have chosen the numerical parameters 
such that L= II = 4D = 1 and we have set d = 0. In part 2 we consider in more detail tlie dependence 
of the deformation properties on the fault parameters. In order to appreciate the magnitude of 
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dcronnatious it will ocoassionally be couveuietit to set 2L= 100 kmaiul to set Uq = I m although 
scaling to other values of L and u^j is straight forward. We will also have occasion to compare the 
present results with those obtained for an Infinitely long screw dislocation (Cohen, 1979), Here we 
again take H = 4D = 50 km but now, of course, L = <», The illustrative discussion is for model DILE. 
For the siike of brevity we will not explicitly discuss tlie results for model DILV although the rele- 
vant contours of displacements and strains will be presented. We warn the reader that the results for 
our two models arc often very different as the postscismic relaxation proceeds. We liuvc chosen to 
concentrate on model DILE since at least some materials appear to remain dilatationally elastic even 
when deviate rically viscoelastic. 

Contours of equal displacements for DILE arc shown in Figures 3 and 4. Strain contours 
are presented on Figures S tlirough 7 and tl\e rotations on Figure 8. The corresponding results 
for DILV are shown on Figures 9 through 13. (Tlie rotations are the same in tlic two models.) 

A useful way to begin tlie interpretation of these results is to consider the variation of 
displacement, ii, and shear strain, ei 2 ) with distance from the fault along the y axis. Along this 
axis the effects of the fault edge are least noticeable and in fact v, ej i , and £32 which depend 
on finite fault sizes for their existence vanisli wlicn x = 0. Even on the y axis, liowever, some 
effects of finite fault sizes will be seen at distant points from the fault. Figure 14 shows the 
expected inonotonic decrease in displacement with distance from tlie fault. The lithospheric 
relaxation reduces the displacements witli the maximum ciiange occurring between y = 25 km 
and y = 50 km. The asthcnosplieric relaxation produces an opposite postscismic motion wlilch 
increases the displacement. The net ofl'eet of both relaxations is to produce a maximum forward 
postscismic displacement of about 3 nun near y = 25 kni and a maximum reverse displacement 
of about 1 mm over a wide zone near, say, y = 175 km. A comparison of these results with 
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those for tlic iiiruntc fault can be curried out with the aid of Figure I4b. As expected for fixed 
H and D tlie displacements at distant points from the fault arc larger wlien L is infinite. Furtlier- 
more, tiic displacements induced by astlicnosphcric relaxation become significant relative to the 
coseismic ones at points closer to the fault when the finite fault size is considered. 

The postseismic displacements off the y axis are complex, particularly near the fault edges 
where there is a coiwergcnce of contour lines. One curious feature associated with the litiio- 
sphcric relaxation is a small region of increased displacements beyond the fault edge. 

Returning again to tlie y axis and examining Figure 15 we find that tire shear strain drop is 

about 10~^ at the fault, fulls to zero near 40 kilometers, and tlien becomes positive taking on a 

value of a few parts in 10® year near y=70 km, before decreasing. This reversal in the sign of the 

shear strain is a clear demonstration of the importance of fault finite size. For an infinite fault 

only the term — contributes to the shear strain. Since u decreases monotonicatly with distance 

3v 

from the fault the shear strain is always negative. In the finite Fault model the term must 

ox 

also be considered, it can become the dominate term at certain points not too close to tlie fault 
and produce a positive shear strain. As for the postseismic strain, tlie lithospheric relaxation 
reduces the strain near tiie faidt. At more distant points tiie postseismic strain is positive and at 
still more distant points it becomes negative again. These strain changes arc fractionally small 
near the fault but significant in regions wlicre the coseismic strain vanislies. Tlie astlicnosplicric 
relaxation produces a broad region of strain increase, 

Leaving the y axis tlie behavior of the postseismic strain is complex with regions of both 
positive and negative strain cjiangcs apparent. Lithospheric relaxation produces a small shear 
strain increase along tlic .x axis near the fault and beyond |x| = L. Conversely, there is a strain 
decrease in a region between the x and y axes. The asthcnosphcric relaxation decreases tiie 
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shear strain near tlic fault and along botli coordinate axes but increases it at sonic distant points 
between axes, 

Tlic perpendicular displacements, v, and the normal strains Cjj and Cjj share the common 
feature of being focussed on the edges of the slipped region. Tlic v displacements tend to de- 
crease as the lithosphere relaxes although Figure 4 docs reveal a region of increased displacement. 
Tlic asthenosphcric adjustment increases the displacements. The net postscismic motion is a de- 
creased displacement near the fault and an increased displacement at distant points. 

Tlic normal strains can involve shortening or extension depending on location. The previously 
observed general tendency for lithospheric relaxation to bring contours closer together and for 
astlienospheric relaxation to separate tiieni can be observed for these strains. When the postscismic 
behavior is examined in more detail close to the fault we find that the initial change in Cj j is a 
shortening (in the first quadrant). As the asthcnospJicrc relaxes there is an extension but the net 
motion remains a shortening. The postscismic behavior of c -22 ”^ar the fault is an extension both as 
the lithosphere and asthcnospherc relax. We should re-cinphasize that th results presented here 
arc for the case L = H = 4D. An analysis of the dependence of the time dependent deformations on 
the choice of fault parameters will be included in Part 2. 

SUMMARY 

Tills paper has considered litliosiilieric and astJicnosphcric viscoelasticity as causes of postscismic 
surface deformations. Displacement and strain equations have been derived for a strike-slip fault 
and illustrated for a particular case involving a moderately shallow fault that cuts the surface, a some- 
what deeper asthenosphcrc-lithosphcrc interface, and a finite fault length. The calculations suggest 
that near the fault the postscismic deformations are significant for moderate and large earthquakes 
although tlic deformations are only a small fraction of the coscisniic ones. At various iioints away 



from the fault, the deformations produced by lithosplicric relaxation exceed the coseismic deforma- 
tions; furthermore, over broad distant regions the slower deformations produced by astltcnospheric 
relaxation also exceed the coseismic ones. Tiie sense of direction of these postscismic motions not 
only differs from the sense of the coseismic motions but also varies witlt location and time. 

Accompanying the postscismic changes in strain arc changes in stress. However, as the litho- 
sphere relaxes or softens tlie reduction in its effective moduli causes stress changes that may be very 
much different from that expected from strain changes alone, Tlie time-dependent changes in stress 
will be studied in detail in Part 2. 

We wish to emphasize tliat we do not regard viscoelastic motions as tlie only possible source of 
postscismic deformations. Other factors, such as postscismic fault creep and timc-ilcpcndent rock 
failure, may be Important sources of deformation. 
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Table I 

Symbol Definition for Displacement and Strain Equations 


X, y,z 

L 

D 

d 

H 

P - 
q = 

P+ = 
= 

% = 
q> = 

t 

Si^ = 
s/ = 

Su' = 
Si.2 = 
52.^== 
S2_2 = 
u 

^0 

V 

M] “ 

a = 
b = 
c = 


coordinates of observation point 
coordinates of a point on the fault 
half-fault length = 2 L -L < f < L 
depth of lower boundary of fault 
depth of upper boundary of fault 
depth to lithosphere-asthenospherc boundary 
f + zj 
f “ z 


d < f < D 


for z = 0, p = q 

f + z + 2mH( 
f - z + 2 mH 
^ - z - 2mH 
f + z + 2 mH 
S-x 


for z = 0, p^ = p_ 
for z = 0, 


forz = 0 , Sj =Sj =S 
for z = 0, S4 = Si^ = S'! 
for z = 0, = S2_ = S’j 


t 3 + y2 +q2) 
f2 + y2 ^ p2 j 

t2 +y2 +q^2 
t2 +y2 +q_2 
t^ + y2 + p^2 
t^ + y^ + pj- , 
displacement parallel to fault 
dislocation across fault 
displacement perpendicular to fault 
bulk modulus of lithosphere 
shear modulus of lithosphere 
shear modulus of asthenosphere 
( 3 k^ +/ii)/( 3 ki + 4 /ij) 

(3kj -Mi)/(3ki +Mj) 

3/^j (3kj -2Mi)/(3k^ +n^) 


Tabic II 

Fffcclive Moduli of Lithosphere and Asthenosphere 
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Tabic III 

Coefficient in Displacement and Strain liquations; Equations 14, 15, I7~20; 


•*1 

3k + 1 
3k 4 


3 

3k + 1 


9k^ + 1 5k - 5 
(3k ♦ I)(3k +4”) 


3 

3k + 1 

“3 

3k + 7 
3k 4 


3k + 1 
3k + 4 




9k^ + 1 5k - 5 
(3k + 1 K3k + 4) 


6k + 2 
3k + 4 


3k + 7 
3k + 4 


1 8k^ + 30k - 1 0 
(3k + IK3k+4) 


6 

3k + 4 

^7 = 

12 

h 

IH 

3k + 4 

'’7 

(3k + 1 )(3k + 4) 

3g : 

6k + 2 
3k + 4 

*’8 

I8k^-6k-22 
(3k + IK3k + 4) 

a, : 

6k + 2 
3k + 4 

b. 

6k + 2 
3k + 4 

^10 : 

2 

^to 

6k + 2 

~ jir ^ 

at, : 

7 

b,, 

6k + 2 
3k + 4 


1 

h,2 

I8k^+30k-I0 
(3k + I)(3k +4) 
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Of POOR QUAUTY 


Figure 1 . Model of a Strike-Slip Fault in a Layer Lying Over a Half-Space. 




(A) Rheological F.lemcnt for Standard Viscoelastic Solid. (B) Rhiological 
FlcmL'nt for Maxwell Substance. 
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Figure 6, Normal Strain, 6i|(t), for Model DILE; 







Figure 8. Rotation cu(t) for Models DILE and DILV; 
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Figure 15. (A) Sliear Strdin, 6)20), Versus Distance from the Fault for Model DILF;Ug = 1 m, L = 50 km for the Scales on Top and 
Right Side of Figure. (B) Shear Strain. Ci 2 fO, Versus Distance from the Fault for Infinite Fault; u© = 1 m. L = 50 km for 
the Scales on Top and Right Side of Figure. 





